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We present a systematic study of thermal conductivity, specific heat, electrical resistivity, ther- 
mopower and x-ray diffraction measurements performed on single-crystalline samples of the pseudo- 
quaternary type-I clathrate system SrgGaieSiao-^Ge^, in the full range of < x < 30. All the 
samples show metallic behavior with n-type majority carriers. However, the thermal conductivity 
and specific heat strongly depend on x. Upon increasing x from to 30, the lattice parameter in- 
creases by 3%, from 10.446 to 10.726 A, and the localized vibrational energies of the Sr guest ions in 
the tetrakaidekahedron (dodecahedron) cages decrease from 59 (120) K to 35 (90) K. Furthermore, 
the lattice thermal conductivity at low temperatures is largely suppressed. In fact, a crystalline peak 
found at 15 K for x = gradually decreases and disappears for x > 20, evolving into the anomalous 
glass-like behavior observed for x = 30. It is found that the increase of the free space for the Sr 
guest motion directly correlates with a continuous transition from on-center harmonic vibration 
to off-center anharmonic vibration, with consequent increase in the coupling strength between the 
guest's low-energy modes and the cage's acoustic phonon modes. 

PACS numbers: 72.15.Jf, 72.20.Pa, 82.75.-z 



I. INTRODUCTION 

Semiconducting clathrate compounds are attracting 
considerable attention because of their potential for ther- 
moelectric conversion applications <A The efficiency of a 
thermoelectric material at a given operation temperature 
T can be quantified by the dimensionless figure of merit 
ZT = S 2 T/ p{n e i + kt), where S, p, n e i, kl, are the ther- 
mopower, electrical resistivity, electronic thermal con- 
ductivity and lattice thermal conductivity of the mate- 
rial, respectively. Intermetallic clathrates are compounds 
consisting of polyhedral cages (basically formed by Si, 
Ge and Sn through diamond-like bonding) that are nor- 
mally filled with monovalent or divalent guest cations. 2,3 
Many of them follow the Zintl rule, where the cage 
atoms are partially substituted by acceptor atoms for 
charge compensation between guests and cages. In ad- 
dition to the large S(T) and small K e i(T) often observed 
in Zintl materials, the most pronounced feature of the 
clathrates is their very low lattice thermal conductivity 
kl (of order 1 W/m K at room temperature) Some 
of these compounds even show glasslike temperature- 
dependant thermal conductivity, although they crystal- 
lize in well-defined structures. Therefore, clathrates 
are good candidates to fulfill the phonon-glass electron- 
crystal (PGEC) concept,^ which is a guideline to search 
for high-performance thermoelectric materials with the 
very rare combination of simultaneously low thermal con- 
ductivity and electrical resistivity. The determination of 
which mechanisms are dominant in lowering kl(T) in 
clathrates presents motivation from the performance im- 
provement perspective, as well as from that of further 
understanding the physics behind atoms vibrating in un- 
conventional crystalline lattices. 

Among the several possible structures formed by these 
materials, the type-I clathrate structure adopted by the 



A 8 Gai6Ge3o group (A = Ba, Sr and Eu) has shown par- 
ticularly favorable thermoelectric properties ji^i The 
unit cell of this bec structure (Pm3n, No. 223) consists 
of 46 cage atoms arranged in two dodecahedra and six 
tctrakaidccahcdra, which incorporate two A atoms at the 
2a site and six A atoms at the 6d site, respectively^! 3 . Ini- 
tial investigations of AgGai6Ge3o clearly indicated that 
the kl below room temperature was lowered in direct 
relation to the decrease of the guest atom's ionic ra- 
dius (Ba +2 largest, Eu +2 smallest), whereas the lattice 
parameter of the three compounds remains relatively 
unchanged.^ This suggested that the general lowering in 
thermal conductivity is not so much related to the mass 
of the guest ions, but rather realized by the fact that 
they are loosely bound to an oversized cage, giving rise 
to an anomalous vibration showing low-frequency, non- 
dispersive localized modes (rattling) which couple to the 
heat-carrying acoustic phonon modes of the rigid cage 
structure and scatter them efficient ly£ & 1 1 1 1 1 1 2 

In addition, n{T) for Ba§Gai6Ge3o exhibited a large 
peak at 15 K, being characteristic of a crystal latticed 
By contrast, k(T) for SrgGaigGeao and EusGai6Ge3o 
showed all of the characteristics of a structural glass4££ 
Neutron diffraction measurements^ revealed that the Ba 
atom in the large cage is located essentially at the center 
of the cage (6d site), whereas a substantial probability 
exists for the Sr atom to move off the site center about 
0.3 A, to one of four crystallographically equivalent po- 
sitions (24j or 24fc sites), and Eu atoms move away even 
more, 0.4 A from the 6d site, suggesting that off-center 
rattling may be necessary to produce glass-like ther- 
mal conductivity i 13 ' 14 i 15 ' 16 At lower temperature, nu- 
clear tunneling among the four site o 17 ' 18 i 19 may also play 
the role. The on-center vibrational freedom of the Ba 
ions can be adequately described assuming independent 
harmonic oscillators (Einstein model) ^ i 13 ' 20 i 21 ' 22 but the 
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TABLE I: Starting (flux) composition, crystal composition, lattice parameter a, electrical resistivity p, thermopower S and 
carrier concentration n at room temperature of SrsGaigSiao-ajGea; 





starting composition 


Sr 


Ga 


Si 


Ge 


a (A) 


p280ir(mfi cm) 


S2S0K (A»V/K) 


n (10 20 /cm 3 ) 


x=0 


8 : 38 : 30 : 


8 


13.6 


32.4 





10.446 


0.26 


-13 


46 


5 


8 : 38 : 24 : 6 


8 


13.7 


27.0 


5.3 


10.483 


0.28 


-26 


40 


20 


8 : 38 : 15 : 15 


8 


15.7 


10.0 


20.3 


10.638 


0.68 


-60 


16 


26 


8 : 38 : 6 : 24 


8 


15.9 


4.0 


26.1 


10.703 


1.23 


-133 


10 


30 


8 : 38 : : 30 


8 


15.9 





30.1 


10.726 


1.85 


-200 
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Sr and Eu vibration cannot be satisfactorily modeled 
this way, indicating that anharmonic vibration contribu- 
tions gain significance in these casesJ^ For SrgGai6Ge3o, 
the anomalous specific heat contribution of the Sr ions 
was successfully reproduced by using a soft-potential 
modeli^i The appearance of a plateau in k(T) and a 
broadened maximum in C/T 3 were interpreted as resul- 
tant from the low-energy excitations associated with an- 
harmonic, quasi-localized vibrations of the Sr ions in the 
cage. 

As an added degree of complexity, later investiga- 
tions on several Ba-based clathrates have shown a strong 
dependence of the thermal conductivity on the major- 
ity charge carrier type ; 20 ' 23 such that compounds with 
n-type carriers show low-temperature crystalline peaks 
in k(T) while their p-type counterparts show signifi- 
cantly lower and glasslike n(T) in the same tempera- 
ture interva h 20 ' 22 i 23 ' 24 i 25 These results appear inconsis- 
tent with the idea of the guest ion vibration having the 
relevant role in producing the glass- like behavior. An 
alternate mode l 20 i 24 ' 26 i 27 based on phonon-scattering by 
charge carriers^ was proposed instead, and the question 
of which factors are dominant at low temperatures re- 
mains a currently open debated 

Given this scenario, it seemed instructive to pur- 
sue an investigation of the guest vibrational behaviors 
and their effect on thermal conductivity from the oppo- 
site approach, by maintaining the guest atom and car- 
rier type fixed, while varying the cage environment in 
a controlled and systematic manner. The compound 
SrgGai6- a Si3o+ a (0 < y < 5) has only been character- 
ized structurally so far— i 29 i 30 and was reported to adopt 
the same type-I clathrate structure as A 8 Gai6Ge 30 . The 
lattice parameter is somewhat dependent on y due to 
the size difference between Ga and Si, but remains be- 
tween 10.460 A (y = 0) and 10.408 A (y = 5), there- 
fore about 3% smaller than that of SrgGai6Ge3o.— Fur- 
thermore, partial solid solutions of Si-Ge in a clathrate 
structure have already been realized in polycrystalline 
Ba§Gai6Si30— x 

Ge x i 31 i 32 Based on all the aforementioned, 
the pseudo-quaternary system SrgGaigSiao-zGez ap- 
pears to offer a suitable opportunity for investigating 
the relationships between cage size, guest vibration and 



thermal conductivity. Indeed, we will show that homo- 
geneous single-crystalline samples within the full range 
< x < 30 can be achieved and, upon increase of x and 
expansion of the cage size, we find a continuous evolu- 
tion from on-center harmonic vibration and crystal-like 
kl(T), to off-center anharmonic vibration and glass- like 
kl(T). 



II. EXPERIMENTAL 

Polyhedral single crystals of 1-5 mm in diam. were 
grown from a self-flux method using excess Ga. High 
purity elements were mixed in an atomic ratio of 
Sr:Ga:Si:Ge = 8:38:(30-X):X (X = 0, 6, 15, 24, 30) in 
an argon filled glovebox. The mixture was sealed in an 
evacuated and carbonized quartz tube, soaked at 1180 °C 
for 2-3 h, cooled over 10 h to 1000 C and then slowly 
cooled over 100 h to 800-700° C. At this point the am- 
poules were quickly removed from the furnace and the 
remaining molten Ga flux was separated by centrifuging. 
The composition of the crystals was examined by electron 
probe microanalysis (EPMA) with a wave length disper- 
sive JEOL JXA-8200 system. The results are shown in 
Table I. The compositions of Si and Ge in the crystals 
are somewhat different from the starting compositions, 
and Ga is deficient in the Si rich crystals. Hereafter, 
the Ge content in the crystal x is used to denote the 
SrgGaigSiso-xGea; samples. 

The crystal structures were refined with powder x-ray 
diffraction (XRD). The spectra of all samples shown in 
Fig. [T] were indexed on the basis of the type-I clathrate 
structure. Since both the atomic size and bond length of 
Ge are larger than for Si, the lattice parameter linearly 
increases as x increases from to 30 in the series of sam- 
ples, such that there is a 3% difference between the end 
compounds (see inset). For x = 20 and 26, higher angle 
peaks like [530] and [611] are visibly broadened, possibly 
due to the disorder among Si, Ga and Ge atoms. Elec- 
trical resistivity, thermopower and Hall coefficient were 
measured in home-made systems by a standard DC four- 
probe method, differential method and DC technique, 
respectively in the temperature range from 4 to 300 K. 
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FIG. 1: Powder x-ray diffraction spectra (Cu Ka radiation) 
for the SrsGaieSiao-ajGea; specimens, and refined lattice pa- 
rameters (inset). 



Thermal conductivity experiments were performed using 
a steady-state method in a home-made cryostat. The 
data is reliable up to about 150 K, above which the ef- 
fect of thermal losses by radiation and wire conduction 
require corrections. The specific heat from 0.3 to 300 K 
was measured using a Quantum Design PPMS with a 
thermal-relaxation method. 
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FIG. 2: Temperature dependence of electrical resistivity p 
and thermopower S for SrsGaieSiso-zGea;. 



III. RESULTS AND DISCUSSION 

The values of electrical resistivity (p), thermopower 
(5), and carrier concentration (n) at room temperature 
are listed in Table I. The latter was estimated from the 
Hall coefficient assuming one type of carriers. The de- 
creasing trend in n with increasing x is consistent with 
the systematic increase of p (T = 280 K) from 0.26 to 
1.85 mfl cm. The temperature dependences of p and S 
are shown in Fig. [2] The monotonic decrease of p(T) 
upon cooling is characteristic of heavily doped semicon- 
ductors or low carrier density metals. For all samples, 
S is negative and the absolute value increases monoton- 
ically with the increase of x. These variations of the 
transport properties can be attributed to a systematic 
decrease in n, from 46 x 10 20 /cm 3 to 3.6 x 10 20 /cm 3 . 
This arises from the fact that Si-rich s amp les tend to be 
Ga deficient (see Table I and refs. 2,29 30) while Ge-rich 
samples follow the Zintl rule more strictly, and are there- 
fore less electron-doped than Si-rich samples. 

The specific heats are plotted as C/T 3 vs. T in Fig. [3] 
In this plotting style, the contribution of rattling ions 
appears as a broad peak over an electronic and Debye 
"background" from the stiff cage. Upon cooling, C/T 8 
for all samples initially rises into the broad peak, followed 
by a local minimum, and finally rises again as 1 /T due to 



the contribution from conduction electron and/or tunnel- 
ing of guest ions. With increasing x, the peak height rises 
systematically and its temperature shifts from 10.5 K to 
7 K. This is not related to n, whose contribution becomes 
vanishingly small at temperatures above 4 K. Rather, it 
is already a qualitative demonstration of how the rat- 
tling of the Sr ions increases in direct relation to the 
cage size. Further information about the Sr vibration 
characteristics can be obtained through a more careful 
analysis of the data as follows, using the same method- 
ology we previously developed to analyze the data of C 
of BasGaieGeao^ 2 - 

As a first approximation, the Sr atoms can be consid- 
ered independent Einstein oscillators, and the framework 
composed of (Ga,Si,Ge) 46 cages a stiff Debye solid. Fol- 
lowing this approach, C of SrgGaieSiao-ajGea; is treated 
as a sum of three terms: an electronic contribution C e i, 
a Debye contribution Cd, and an Einstein contribution 
Ce- As we described in detail in ref.[22], first the Sommer- 
feld coefficient 7 and the Debye temperature 0d should 
be evaluated independently and fixed, together with the 
pre-defined dimensionalities and numbers of Einstein os- 
cillators, so that only two fitting parameters are left, 
which characterize the Sr guests vibrational energies. It 
should be recalled that the six Sr( 6d ) ions show a strongly 
anisotropic vibration with greater amplitude within the 
plane parallel to the larger cage's two hexagons.— Because 
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FIG. 3: Temperature dependence of specific heat C for 
SrgGaieSiso-iGe^, presented as C/T vs T. 

the dimensionality plays a role in the Einstein model, 
at least two vibrational modes should be required to 
describe the 6d site alone: in-plane 9 E ^ d ^ and out-of- 
plane 9 E ^ d y respectively^ In addition, a third vibra- 
tional mode E (2a) is required to account for the smaller, 
but still Einstein-type rattling of the two Sr( 2a ) which can 
be assumed isotropic^ In our model, the dimensionalities 
and numbers of oscillators are predefined: N E (2a) =3x2, 
^E(6d) =2x6, = 1x6. We impose the addi- 

tional constraints 9 E{6d) < d^ {M)1 6 E{2 a), E {6d) = ° EL 
and 9 E (Q d j = 0£(2a) = &eh, so the fitting parameters are 
only the lower and higher Einstein temperatures, 9 el 
and 9 E H respectively. 

From C/T vs. T 2 plots (not shown) the obtained val- 
ues of 7 are between 11 and 24 mJ/mol K 2 , and the 9o 
decreases from 370 K to 200 K on going from x = to 
30. However, as we discussed previously,— when Sr^d) 
anharmonic vibration becomes relevant it interferes with 
this evaluation even at the lowest temperatures. In fact, 
this last value of 200 K for 9jj of Sr 8 Gai 6 Ge3o is an arti- 
fact, much smaller than 312 K estimated from atomic dis- 
placement parameters^ Thus, we use the value of 288 K 
obtained for BagGai6Ge3o as a better representation of 
the Ga-Ge cages' Debye temperature for SrgGai6Ge3o. 
For x = 5, 20, 26, 0d is estimated by linear interpolation 
between the values for x = and x = 30. 

Fig. H shows the fits to the data of C/T 3 for x = 0, 
20 and 30. For x — 0, an excellent fit is obtained with 
9 el — 59 K and 9eh = 120 K, indicating that the as- 
sumption of independent, harmonic oscillators contained 
within the Einstein model applies very well for the Sr vi- 
bration in SrgGaigSiao. By x = 20 it is clear that the fit 
with the Einstein model is no longer a good representa- 
tion for the data, and at x = 30 the model has become 
completely inadequate to reproduce the behavior. This 
is a clear indication of the effects caused by increasingly 
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FIG. 4: (Color online) Einstein model fits (green solid lines) 
of the C/T data (symbols). The dotted, dashed and dashed- 
dotted lines are the calculated Ce, Cd, and C e i, respectively 
(see text). 



anharmonic vibration of the Sr^) guest ions, as their 
potential well broadens in direct proportion to the free 
space available for movement J£ In such soft 
potential model gives a good fit to C(T)£L To our knowl- 
edge, however, there are no theoretical models available 
that can describe the C(T) in this continuous transition 
from harmonic to anharmonic vibration observed in our 
sample series. For x = 30, the value o^9el = 35 K deter- 
mined mainly by the peak position is in good agreement 
with the result of 38 K observed by Raman scattering.— 
The Einstein temperatures 9 el and 9 eh for all samples 
are listed in Table II. 

Let us now see how all of these systematic changes 
affect the heat transport, which is our main purpose. 
The total thermal conductivity k of SrsGaigSiso-^Gea, 
is plotted as a function of temperature up to 150 K in 
Fig. O We can see two types of systematic evolutions in 
the data with increasing x: (i) n at higher temperature 
decreases by a factor of 3, and (ii) a low temperature 
crystalline peak for x — is gradually but strongly sup- 
pressed, resulting in glasslike behavior for x — 30. The 
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TABLE II: Parameters used for the solid line curves in Fig. [4] 
and Fig. HJ] 



Parameter Unit x=0 5 20 26 30 



0EL 
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59 
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40.5 
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6eh 
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120 


115 


100 


94 


90 


Ci 


10 30 /(m s 2 K 2 ) 


0.3 


0.5 


1.3 


2.1 


4.2 


c 2 


10 30 /(m s 2 K 2 ) 


0.1 


0.16 


0.4 


0.7 


1.4 


A/B 


10 5 /(m K) 
l/(m K 4 ) 


0.65 


1.1 
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4 
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D 
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FIG. 5: Temperature dependence of total thermal conductiv- 
ity k. The inset shows the estimated electronic part K e i. 
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At lower temperatures, the plot shows a clearer picture 
of the peak suppression. We will analyze the data using 
the same approach described in detail in rcf. For lack 
of a single model capable of describing all sample behav- 
iors in C(T), we use the results of the Einstein model fit- 
tings despite their poor quality in Ge-rich samples. This 
is justifiable in our comparative analysis because kl(T), 
given in the semiclassical theory by 

*L = \j ° du[C L {u,T)vl], (1) 

is not limited by the lattice specific heat Cl(oj, T), nor 
the average sound velocity v, but rather by the very low 
phonon mean free path I, which is averaged over all major 
contributing scattering mechanisms. In the TRR model/ 
I is written as 



FIG. 6: Temperature dependence of lattice thermal conduc- 
tivity kl- The solid lines are fits using the TRR model de- 
scribed in the text. 



first effect can be directly attributed to the electronic 
contribution K e i(T), which is estimated from the elec- 
trical resistivity p(T) using the Wiedemann-Franz law, 
Kel (T) = (n 2 k 2 B /3e 2 )T/p(T). As shown in the inset of 
Fig. K e i (T) have a simple, slightly sublincar behavior 
with slopes directly related to the carrier concentration, 
and therefore no relevant effect on the low temperature 
peak. 

By subtracting K e i(T) from k(T), we can estimate 
the lattice contributions kl(T), which are shown in 
Fig. [SJ The values are now much closer at tempera- 
tures above 100 K, and an interesting feature reveals itself 
in this range: the heat conduction level of the pseudo- 
quaternary (intermediate) samples is lowered with re- 
spect to the ternary (end) samples. This is most likely 
the effect of extra phonon scattering on Ga/Si/Ge site 
disorder on the cage. 



' — i^TS + ^res + ^ Ray) + Imin (2) 

which includes three mechanisms: tunneling between 
localized guest sites, resonance scattering from guest ion 
rattling and Rayleigh scattering from impurities, imper- 
fections and mass difference. 

The best fits to the data are shown as solid lines in 
Fig. [S] and the parameter values are summarized in Table 
II. The most relevant results are the increase by one order 
of magnitude in both the resonant scattering level Ci 
and the TS scattering level as x increases. The latter 
can be expressed by the ratio A/B = n(hv) 2 /fcs, which 
in glasses is essentially a measure of subset density of 
tunneling states h that are able to strongly couple to the 
phonons and effectively scatter them^ For x — 5, 20, 26, 
Si disorder among Ga and Ge makes Rayleigh scattering 
level larger than that for x — and 30. 

Thus, the combined results of lattice Cl(T) and kl(T) 
can be clearly attributed to a systematic evolution of the 
Sr/g^N guests rattling level, with a decrease in character- 
istic energy and increase in anharmonicity, both arising 
from the fact that the lattice expansion increases the free 



6 




1.25 1.30 1.35 

guest free space (A) 



1.40 



FIG. 7: Dependence of the lower Einstein temperature 9 el 
and lattice thermal conductivity at 15 K on the guest free 
space. 



space for guest excursion and deforms the restoration po- 
tential. We may conclude that this is the main cause of 
the peak suppression in kl, leading to glasslike behavior. 
Fig. [7] summarizes these results by plotting kl at 15 K 
and 6 el in terms of the guest free space, evaluated semi- 
quantitatively by subtracting the Sr ionic radius (1.35 A) 
and Ga covalent radius (1.26 A) from the tetrakaideca- 
hedron's radius in each sample. 



clathrate structure's cage size without changing the 
guest ion or the charge carrier type. The lattice param- 
eter expands by up to 3% with increasing Ge content, 
and as a consequence the free space for guest excursion 
increases in the cage. The characteristic energy of the 
localized Sr( 6d ) vibrations decreases from 59 K to 35 K, 
and the Sr( 6d ) behavior clearly evolves from vibrating 
in an on-center harmonic potential in SrsGaigSiao to 
a broadened potential in SrgGai6Ge3o, which allows 
off-center and anharmonic vibration. This leads to 
an increase in the effective density of tunneling states 
and a strong enhancement of the coupling between 
the Sr(6d) vibration and the cage acoustic phonons, 
shortening the latter's mean free path. As a result, the 
low temperature (1-20 K) lattice thermal conductivity 
is suppressed in such a way that the crystalline-like 
peak found for SrgGai6Si3o evolves into the well-known 
glasslike behavior of SrgGai6Ge3o. 

Our results leave little doubt that the described mech- 
anism is the dominant one in producing the anoma- 
lous thermal conductivity behaviors observed in these 
clathrates at low-intermediate temperatures. However, 
other factors certainly contribute, with variable relevance 
depending on the particular system or temperature range 
under study. 
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